11" 


JULY  1963 


LETHAL  AREA  DESCRIPTION 


K.  A.  Myers 


RDT  &  E  Project  No.  1M023201A098 

BALLISTIC  RESEARCH  LABORATORIES 


i 

I 

mmseaBSi 

newi 

«fi»eKaK8s»aQffiMB8S»aMs^ 

MMBBnaPw 

BffiscsmBa 

tKfiiaSSWM 

»9Bffi0SBQ80l 

fi 

s 

laaa 

aSBKSSggBSgSSBfi^^ 

!WS9B«H6H 

nssasssss 

RBSHBSa 

bSsSbSSw 

ABERDEEN  PROVING  GROUND,  MARYLAND 


^  0  o  M  Si  ^  ^ 


TfflS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


REPRODUCED  FROM 
BEST  AVAILABLE  COPY 


*  # 

BALLISTIC  RESEARCH  LABORATORIES 
TECHNICAL  NOTE  NO.  I510 
JULY  1963 


LETHAL  AREA  DESCRIPTION 


K.  A.  N(yers 


Weapon  Systems  Laboratory 


RDT  eo  E  Project  No.  IMO2320IA098 


ABERDEEN  PROVING  GROUND,  MARYLAND 


•  • 


BALLISTIC  RESEARCH  LABORATORIES 
TECHNICAL  NOTE  NO.  I5IO 

KAMyers/  jdk. 

Aberdeen  Proving  Ground,  Md. 
July  1963 


LETHAL  AREA  DESCRIPTION 

ABSTRACT 

A  detailed  description  of  the  lethal  area  problem  is  given, including 
discussion  of  the  procediire  presently  used  at  BRL  in  the  computations  of 
lethal  areas . 
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INTRODUCTION 


This  technical  note  was  originally  prepared  in  1957  as  an  internal 
working  paper.  The  intent  was  to  present  a  detailed  description  of  the  lethal 
area  concept  at  BRL  and  the  method  used  to  compute  lethal  areas  on  high  speed 
digital  computers.  In  view  of  the  wide  interest  in  this  area  over  the  years 
and  as  a  consequence  of  several  recent  requests  for  copies^  this  technical 
note  has  been  prepared. 

The  initial  portion  of  this  paper  consists  of  an  up-dating  of  the  lethal 

(l)* 

area  concepts  discussed  in  BRL  Report  No.  oOO  which  was  necessitated  by  the 
introduction  of  more  refined  casualty  criteria  and  presented  area  functions. 

The  remainder  of  the  paper  is  devoted  to  a  discussion  and  description 
of  the  techniques  and  procediires  presently  used  in  lethal  area  evaluations. 
Although,  the  possible  extensions  of  the  problem  are  not  covered  specifically, 
it  is  felt  that  sufficient  background  is  given  to  allow  these  special  cases 
to  be  handled  by  the  reader  as  necessary. 

LETHAL  AREA  CONCEPT  AND  INPUT  REQUIRED 

Lethal  area  can  be  expressed  mathematically  as  follows:  Define  d(x,y) 
as  the  density  of  targets  in  an  element  of  area  about  the  point  (x,y)  and 
denote  the  probability  that  a  target  in  that  element  of  area  will  be  incapaci¬ 
tated  as  P  (x,y).  We  may  then  write  the  equation  for  the  expected  number  of 

IV 

casualties,  E^,  as  follows: 

-h»  -Ko 

E^,  =  /  /  a(x,y)  Pj^(x,y)  dxdy. 

-00  —00 

Now  if  we  assume  the  targets  are  uniformly  distributed  over  the  ground  plane, 
CT(x,y)  can  be  represented  simply  as  a  constant,  a,  and  we  can  write: 

E  +00  +<» 

■^  =  /  I  dxdy. 

—00  -00 


* 

The  parenthetical  superscripts  refer  to  reference  numbers. 
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The  quantity  —  has  the  dimensions  of  area  and  is  called  the  lethal 
0 

area.  Some  operations  analysts  often  refer  to  this  quantity  as  the 

mean  area  of  effectiveness  (M.A.E.).  It  is  obvious  from  an  examination  of 

this  expression,  that  the  lethal  area  when  multiplied  by  the  density  of  targets, 

will  yield  the  expected  number  of  casualties.  A  further  examination  of  the 

above  expression  indicates  that  lethal  area  is  really  just  a  weighted  area; 

the  weight  for  each  element  of  area  being  determined  by  the  probability  of 

incapacitation  function  P  (x,y). 

K 

Let  us  now  look  at  some  of  the  geometry  associated  with  the  problem. 

(See  Figirre  l). 


Figure  1 
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Assvune  that  the  shell  is  approaching  the  ground  at  an  angle  of  fall  oo  and 
velocity  The  shell  bursts  at  height  h.  The  problem  is  to  find  the 

probability  that  a  target  at  the  point  {x,y)  is  "killed"  by  the  shell.  We 
see  that  the  target  will  be  located  in  a  small  fragment  spray  leaving  the 
shell  about  an  angle  0  measured  from  the  nose  of  the  shell.  In  order  to 
determine  the  probability  that  the  target  will  be  a  casualty,  a  considerable 
amount  of  information  will  be  required.  Specifically,  the  type  of  information 
required  may  be  listed  as  follows: 

(1)  Velocity  fall  off  law. 

(2)  Casualty  criterion. 

(3)  Target  presented  area. 

(4)  Fragment  mass  distribution. 

(5)  Fragment  density. 

(6)  Initial  fragment  velocity. 


The  manner  in  which  the  velocity  falls  off  with  range,  r,  is  usually  approxi¬ 
mated  in  the  following  manner: 


C,  p  A 


f  r 


V  e 
o 


m 


where  is  the  initial  fragment  velocity^  is  the  drag  coefficient,  p  is 
the  air  density,  A^  is  the  average  presented  area  of  the  fragment,  and  m  is 
the  weight  of  the  fragment.  In  this  connection,  and  A^  are  generally 
determined  experimentally,  vinless  of  course  the  shape  of  the  fragments  are 
such  that  A^  may  be  estimated  theoretically.  In  any  event,  once  these 
constants  have  been  determined,  we  are  in  a  position  to  determine  the  velocity 
of  that  specific  fragment  as  a  function  of  distance  traveled. 

During  the  past  few  years  a  concentrated  effort  has  been  directed  towards 
obtaining  a  realistic  casualty  criterion.  To  this  end,  the  Biophysics  Laboratory 
at  Edgewood  Arsenal  has  developed  information  pertaining  to  the  woimding  power 
of  various  size  fragments  striking  the  body  at  various  velocities.  These  data 
were  analyzed  by  personnel  of  the  Terminal  Ballistics  Laboratory  here  at  BRL. 

The  resulting  casualty  criteria  originally  were  reported  in  BRL  Report  No.'996(^^ 


t 


A  subsequent  medical  reassessment  of  the  basic  wound  data  resulted  in  a  modi- 

(3) 

fication  of  these  criteria  as  described  in  BRL  Technical  Note  No.  1297. 
Essentially,  these  casualty  criteria. are  a  series  of  functions  for  various 
tactical  situations  which  describe  the  probability  that  a  fragment  hitting  a 
target  will  create  a  casualty. 

With  respect  to  the  target  presented  area  A^,  various  "cover  fimctions", 
or  more  properly,  presented  area  functions,  have  been  in  existence  for  some 
time  which  describe  the  presented  area  of  the  target  as  a  function  of  the 
position  of  the  target  relative  to  the  burst.  Some  of  the  more  prominent 
cover  functions  which  have  been  used  are:  (l)  a  variety  of  cut-off  angles 
which  were  intended  to  simulate  foxhole  cover;  (2)  the  ORO  cover  functions 
which  express  the  presented  areas  as  a  function  of  the  aspect  angle;  and 
(3)  most  recently,  some  new  BRL  cover  functions  which  have  been  developed  to 
consider  typical  terrain  features.  The  latter  are  described  in  BRL  Memorandum 
Report  No.  1203!^^ 

The  fragmentation  characteristics  of  the  shell  can  adequately  be  described 
in  terms  of  initial  fragment  velocity,  fragment  density  and  fragment  mass  dis¬ 
tribution.  In  most  cases,  these  data  will  vary  as  a  function  of  0,  the  angle 
measured  from  the  nose  of  the  shell.  Typical  plots  of  initial  fragment  ve¬ 
locity  V^,  and  fragment  density  are  shown  in  Figiore  2. 


Figure  2 
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Fragment  mass  distributions ^  as  the  name  implies,  are  used  to  describe 
the  distribution  of  fragment  sizes  within  a  spray.  These  mass  distributions 
are  also  a  function  of  0.  In  general,  these  fragmentation  characteristics 
are  determined  from  static  detonation  tests  where  fragments  are  collected  in 
recovery  boxes  filled  with  celotex  sheets.  The  fragments  are  then  extracted 
from  the  celotex,  counted  and  weighed.  Fragment  velocity  measurements  are 
made  with  the  aid  of  high  speed  photography.  A  description  of  the  techniques 
employed  in  the  fragmentation  tests  conducted  for  the  BRL  by  Development  & 
Proof  Services  is  given  in  Reference  5* 


How  then  are  these  basic  data  combined  mathematically  to  yield  Pj^(x,y)? 

To  illustrate,  let  us  suppose  that  a  small  spray  about  the  angle  0  is  comprised 
of  only  three  fragment  sizes j  namely  m^  grains,  m^  grains  and  m^  grains,  and 
further,  that  q^  is  the  fraction  of  the  fragments  in  the  spray  weighing  m^ 
grains,  is  the  fraction  of  the  fragments  in  the  spray  weighing  m^  grains  and 
q^  is  the  fraction  of  the  fragments  in  the  spray  weighing  m^  grains.  Since  we 
have  ass\mied  that  the  fragments  in  this  spray  all  leave  the  shell  with  the  same 
velocity,  we  can  employ  the  velocity  fall-off  law  described  previously  and  thus 
associated  with  our  three  fragment  sizes  we  will  have  three  striking  velocities 
Vi,  Vg  and  V^.  Using  the  casualty  criteria  we  can  then  determine  the  proba¬ 
bility  that  the  target  will  be  incapacitated  if  hit  by  a  fragment  weighing  m^^ 
grains  and  moving  at  a  velocity  V^.  We  denote  this  probability  by  P^  ; 

similarly,  we  have  P  and  P  .  If  we  then  describe  the  density  of  all  the 

ixK/-.  nK.^ 

J/ 

fragments  of  the  spray  to  be  fragments  per  steradian,  the  density  of  m^ 


fragments  will  be  qj^q;  fragments  q^q  and  m^  fragments  q^q.  Now,  employing 

A  ,  the  presented  area  of  the  target,  we  should  expect  — 5 —  of  the  m  grain 

^  d. 


fragments ;  — ^ —  of  the  fragments  and  —  of  the  m_  fragments  to  hit 

r  r  ^ 

the  target.  The  probability,  therefore,  that  the  target  will  not  be  "killed"  by 
the  m  fragments  may  be  represented  as  follows : 

q3_Ti  At 


-  ■’hi,) 
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Similarly,  the  probability  that  the  target  would  not  be  ''killed"  by  the  m^ 
fragments  would  be  as  follows: 


and  for  the  fragments  as  follows: 


Then  the  probability  that  the  target  will  not  become  a  casualty  by  either  of 
the  m^^,  nig  or  m^  fragments  can  be  represented  by 


-  X’ 


*  (1  -  Phkj) 


*  <1  -  ) 
P 


For  most  practical  problems,  with  natural  fragmenting  shell,  there  will  be 
many  more  mass  groups  than  I  have  discussed  here,  and  it  would  therefore  be 
advantageous  to  use  the  Poisson  approximation  to  the  binomial  distribution 
Just  discussed,  l.e.: 


^i^At 


Thus,  the  expression  for  the  probability  that  none  of  the  fragments  in  "n" 
mass  groups  would  create  a  casualty  can  be  written: 


-  T)  A 


^  'll  ^hk 

i=l  ^  “‘l 
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We  are  now  in  a 
that  the  target 


position  to  write  the  expression  for 
is  killed,  as  follows: 


the  probability¬ 


'll  ^hk. 


There  has  from  time  to  time  been  considerable  discussion  regarding  this 
"approximation"  of  the  probability  of  incapacitation,  inferring  that  the 
Binomial  approach  represents  a  "true" solution.  Actually  a  more  correct 
application  of  the  binomial  distribution  could  lead  to  an  expression  of  the 
form: 


n 

1  -  Jt  (1  - 

i=l 


t  hk> 


k  r 


where  N  represents  the  total  niimber  of  fragments  in  the  spray  and 
0,  being  the  nimiber  of  steradlans  subtended  by  the  spray. 


This  expression  for  the  probability  of  incapacitation  can  again  be 
approximated  by: 


Even  so,  the  Binomial  itself  must  be  considered  only  as-lapproximation  since 
fragments  projected  from  the  same  shell  cannot  really  be  considered  as 

At  n 

"independent  events".  Further,  if  one  considers  that  — ^  E 

r  i=l  i 

is  the  number  of  incapacitating  fragments  expected  to  strike  the  target,  an 
estimate  of  the  probability  that  the  target  is  hit  by  at  least  one  such 
fragment,  i.e.,  the  probability  that  the  target  is  incapacitated,  can  be 
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estimated  directly  by  the  Poisson  distribution  as: 

C 

O 


^  A  n 

^  '^i  ^hk 

r  i=l  i 


T)  n 

"T“  ^i  V 

r  i=l  i 


^  ^^i  ^hk. 

r  1=1  1 


=  1  -  e 


To  the  author*  s  knowledge,  no  extensive  comparisons  of  lethal  areas  computed 
using  the  Binomial  and  Poisson  forms  for  the  probability  of  incapacitation 
have  been  made.  However,  in  some  isolated  cases  where  such  comparisons  have 
been  made,  the  difference  did  not  appear  significant,  particularly  in  light 
of  other  assumptions  inherent  in  the  lethal  area  concept. 

LETHAL  AREA  COMPUTATION 

Having  obtained  the  input  data  described  in  the  preceding  sections,  one 
is  now  in  a  position  to  proceed  to  the  actual  computation  of  a  lethal  area. 

The  following  sections  describe  the  method  presently  being  used  at  BRL. 

By  a  tremsformation  to  the  coordinate  system  the  lethal  area 

integral  becomes: 

A^  =  /  /  r  Pjj(r,4i)  dr  d®  , 

r  $ 

where  r  is  the  distance  from  the  burst  to  a  point  on  the  ground  and  $  is  an 
angle  in  the  ground  plane  measiured  from  the  projection  of  the  shell  trajectory 
into  that  plane  to  the  line  joining  the  origin  to  the  point  in  the  ground  plane. 

The  minlmvun  distance  to  the  ground  plane  is  h,  the  height  of  burst,  and 
although  theoretically  the  upper  limit  on  r  would  be  infinite  in  practice  we 
shall  use  the  value  r  „  such  that  for  all  r  >  r  ^  the  value  of  the  integrand 


would  be  zero.  Taking  advantage  of  the  symmetry  in  $  the  integral  becomes: 


r  p  Jt 

A,  =  2  /  r  P,  (r,<t)  dr  d®  . 

h  0  ^ 

Or 

r  ,  rt  r  p  Jt 

r  P  (r,$)drdii  +  J  f  r  P  (r,$)drd4i  +  r  f  ^  P,  (r,o)drd$ 

^  h+1  0  ^  r  T  0  ^ 

ml  .  _ 

If  in  the  second  integral  we  make  the  substitution  d  (Ln  r)  =  ^  dr  we  have: 

Ln  r  Jt  p 

Jr  P,  (r,$)  cLi)d(Ln  r). 

Ln(h+1)  0 

It  will  be  noted  that  if  one  integrates  numerically  with  respect  to  Ln  r, 
taking  equal  step  sizes  in  Ln  r,  the  effect  will  be  to  concentrate  the 
evaluations  of  the  integrand  in  close  to  the  point  below  the  burst  where  one 
might  expect  the  integrand  to  change  rapidly,  and  to  spread  out  the  interval 
in  r  between  successive  evaluations  of  the  integrand  at  the  longer  distances. 
It  is  also  apparent  why  the  integration  from  h  to  r^.,  was  separated  into  two 
parts.  Without  the  separation  the  lower  limit  of  the  integral  for  h  =  0, 
would  have  become  Ln  (O).  It  should  also  be  noted  that  the  separation  of  the 
integration  at  h+1  is  purely  arbitrary.  In  practice,  the  first  of  the  three 
integrals  is  approximated  by: 

Jt 

(h+l/2)  /  P  (h+1,  <t)d<l)  . 

0  ^ 

In  evaluating  the  third  integral,  the  integration  in  r  is  accomplished  in 
equal  steps  of  r. 


=  2 


h+1  Jt 

/  I 

h  0 


•  • 


The  lethal  area  integral  as  evaluated  by  the  machine  is  then: 


where  as  indicated  the  integration  over  ®  is  considered  as  part  of  the 
integrand  with  respect  to  r. 

Generally  the  fragmentation  data  are  obtained  for  shell  detonated  stati¬ 
cally.  It  is  necessary  therefore,  to  adjust  the  data  to  simulate  the  fragmen¬ 
tation  characteristics  of  a  shell  bursting  in  flight. 

The  fragmentation  data  obtained  from  tests  consist  of  an  average  initial 
fragment  velocity,  average  fragment  density  and  a  mass  distribution  for  each 
of  several  Intervals  of  the  angle  9  between  0°  and  l80°.  Continuous  functions 
of  initial  fragment  velocity  and  fragment  density  vs.  9  are  obtained  by  plotting 
these  quantities  at  the  midpoints  of  their  respective  intervals  and  connecting 
the  points  by  straight  line  segments.  Mass  distributions  are  considered  to  be 
constant  over  the  interval  which  they  represent. 

Figure  3  indicates  the  scheme  for  adjusting  the  static  fragmentation  data 
to  dynamic  data. 


Figure  5 
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Here  a  fragment,  which  would  have  left  the  shell  at  a  velocity  V  ajid  an 

'  s 

angle  0  had  the  shell  been  at  rest, would  be  thrown  forward  by  the  shell's 
s 

remaining  velocity  to  an  angle  6^^,  the  resultant  initial  velocity  being  Vj^. 
The  following  relations  can  be  used  to  determine  and  0^^. 

1)  V  sin  0  =  sin  0^ 

s  s  D  D 

2)  Cos  =  Vj,  Cos 


In  addition  to  the  change  in  fragment  direction  and  initial  velocity  it 

will  be  noted  that  if  a  small  fragment  spray  of  width  2^  ,  centered  at  0 

s  s 

is  thrown  forward  to  a  spray  of  width  centered  at  0^,  the  density  of 

fragments  in  the  spray  will  change  from  a  to  Cy,.  Since  the  number  of  fragments 
contained  in  the  spray,  N,  will  remain  constant  we  have  the  relations: 


a  = 


N 


N 


Ititsin0  slru^ 
s  s 


sin0j^ 


or 


sln0  siru^fi 
_  s _ s 

D  “  °s  sin0jj  sln^j^ 


then 


Lim  fa  sin0  I  sin  /?i0  I 

IS  s'  si 


1  \  sin0 

d0 

M-  CT  ® 

s 

•1  I  s  sin0^ 

d0^ 

n  j  E 

D 

d0 

g 

In  order  to  compute  we  consider  the  equation: 


D 


V  Cos  0  sin0^  +  V,  sin0T,  =  V  sin0  Cos  0^  . 

s  s  Dh  Ds  s  D 


Now  since  our  function  of  vs.  0  is  made  up  of  a  number  of  straight  line 
segments  we  may  express  V  by: 

5 

V  (O  )  =  A.  e  +  B.  .  0  in  radians 

s'  s'  is  i  '  s 


where  the  constants  A.  and  B.  depend  upon  the  interval  in  which  0  is 

X  X  S 

contained. 


Making  this  substitution  for  V  and  carrying  out  the  differentiation 

s 

we  find: 


To  obtain  the  dynamic  fragmentation  characteristics  needed  for  the  lethal  area 

code  a  supplementary  code  has  been  devised  which  takes  as  input,  a  table  of 

initial  fragment  velocity,  density  and  a  count  number  tabulated  for  each  5° 

from  0  =  0  to  0  =  l80°.  The  count  ntmiber  is  merely  a  device  used  to  identify 

the  mass  distribution  associated  with  a  particular  interval  in  0^.  This  is 

accomplished  by  choosing  count  numbers  so  that  all  count  nimbers  corresponding 

to  angles  of  0  from  a  particular  interval,  when  rounded  off  to  the  nearest 
s 

integer  give  the  Integer  corresponding  to  the  identification  n\miber  of  the 
mass  distribution  associated  with  that  particular  interval. 

The  purpose  of  the  code  is  to  produce  a  table  of  dynamic  initial  fragment 
velocity,  dynamic  density  and  covint  number  tabulated  for  each  5°  of  0j^  from 
0°  to  180°. 

Essentially  what  the  code  does  is  to  determine  what  fragment  angle,  0 

s 

in  the  static  test  would  be  vectored  forward  to  each  angle  0^  as  the  result  of 

a  shell  remaining  velocity  V  .  Having  found  the  angle  0  for  each  0  the 

^  s  u 

machine  interpolates  in  the  input  table  to  determine  the  corresponding  initial 

fragment  velocity,  density  and  count  number.  The  interpolated  count  number 

then  becomes  the  count  number  corresponding  to  the  angle  0^^.  In  this  way  the 

mass  distribution  is  shifted  forward  with  the  0  interval  to  which  it  corresponds. 

Then  the  interpolated  values  of  0  ,  V  and  a  are  used  to  compute  and  (j„. 

s  s  s  D  D 

The  resulting  table  is  used  as  input  for  the  lethal  area  code. 


The  mass  distributions,  each  identified  by  an  integer,  contain  a  series 
of  fragment  weights,  m^,  representative  of  the  i  fragment  weight  group, 
and  a  series  of  normalized  percentages,  representing  the  fraction  of  the 
total  number  of  fragments  in  each  of  the  fragment  weight  Intervals. 

In  addition  to  the  fragmentation  characteristics,  a  table  relating  target 
presented  area  and  aspect  angle,  \lr,  to  the  distance  r  is  needed  for  each 
height  at  which  lethal  area  is  to  be  computed.  The  relation  between,  aspect 
angle,  h,  and  r  is  shown  in  Figure  4. 


Figure  4 

Also  shown  in  Figure  4  are  the  quantities  h,  and  rj  where  h.  is  the 
height  above  the  ground  at  which  it  is  desired  to  compute  the  probability  of 
Incapacitation  and  r'  is  the  distance  which  a  fragment  would  have  to  travel 
in  order  to  hit  that  point.  It  is  apparent  that  by  specifying  the  r,  h,  and 
•(tf,  r'  and  h^  are  fixed.  This  technique  is  used  when  it  is  desired  to  consider 
targets  at  different  ranges  as  being  located  at  different  heights  above  the 
mean  ground  plane  or  where,  as  might  be  the  case  with  a  standing  target,  it 
is  desired  to  base  the  probability  of  incapacitation  upon  the  fragment  spray 
striking  the  center  of  the  man. 


Figure  5  shows  a  sample  input  information  sheet  for  the  lethal  area  code. 
The  following  paragraph  will  describe  the  various  quantities  needed. 


NUMBER  -  INPUT  DATA 


1.  h  = 


R  T  = 
ml 


ft. 

ft. 


= 


Log  Int. 


2.  -K  = 


sin  CD 


M  = 


cos  OD  = 


_ft. 

ft. 


2nd  Int. 


Intervals  cd  = 


-a  = 


max 
J  = 


ft. 


Intervals 


Log  Int. 


Radians 


5.  b  = 


n  = 


1/2 


Radians 


1/5 


It. 


^  = 


Radians  $  =  5.1^159 

max 


ft/ sec  Run  ID 


-.01 


5.  -20 


cutoff  velocity- 


max 


1  X  10 


-8 


No.  of  q  and  m 
_ ^ft/sec 


6.  Vector  Run 

7.  q  and  m _ 

8.  Cover 


FIGURE  5 


1.  h 


initial  height  at  which  lethal  area  is  to  he  computed 

Increment  by  which  burst  height  is  to  be  increased  if  lethal  areas 
are  to  be  rim  for  several  heights 


max 


R 


ml 


^m2 


maximum  burst  height  for  which  lethal  area  is  to  be  computed 
r^i  (see  section  describing  break  down  of  integral) 

r^P  (see  section  describing  break  down  of  integral) 

number  of  intervals  to  be  used  in  integration  from  Ln(h+l)  to 

Ln  r  T  MUST  BE  EVEN  INTEGER. 

ml  . .  '  " 


2.  -K 

M 


from  velocity  falloff  law  =  V^.e 


Kr 

m 


number  of  intervals  to  be  used  in  integration  from  r  ,  to  r  ^ 

ml  m2 

MUST  BE  EVEN  INTEGER. 


sin  CD  CD  =  shell  angle  of  fall 

Cos  CD  ,  CD  =  shell  angle  of  fall 

CD  CD  =  shell  angle  of  fall  MUST  BE  IN  RADIANS, 

-a  casualty  criterion  constant 


3*  b  casualty  criterion  constant 

n  casualty  criterion  constant 

angular  Interval  at  which  fragment  density  and  initial  velocity 
and  "count"  numbers  are  stored  in  machine 


it.  step  size  to  be  used  in  integration  over  $ 

NOTE:  NUMBER  OF  INTERVALS  MUST  BE  EVEN. 

Vj^  shell  remaining  velocity 

Run  ID  number  to  be  used  for  identification  purposes 
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5.  C  number  of  mass  distributions  +0.5  • 

^  max 

Cut  off  velocity  it  is  generally  assumed  that  there  is  some  minimum 

velocity  whicli  a  fragment  must  have^  regardless  of  its 
size,  in  order  that  it  has  any  chance  of  producing 
casualty. 

identification  number  for  table  of  6  ,  V  ,  o  ,  and  count 

s'  s'  s' 

number. 

identification  number  for  whole  group  of  mass  distributions. 
Not  to  be  confused  with  identifying  integers  for  the  indi¬ 
vidual  mass  distributions . 

where  tables  of  presented  area  and  aspect  angle  vs.  r  for  the 
bxirst  heights  of  interest  have  previously  been  given  to  the 
Computing  Laboratory.  An  identifying  symbol  is  used  here. 

At  each  point  (r^,  of  the  integration  grid  the  probability  of 
incapacitation  is  computed  as  follows ; 

The  angle  i|r  and  the  presented  area  A  are  determined  for  the  value'  of  r. 

^  i 

from  the  table  corresponding  to  the  burst  height  being  considered.  The  angle 
©p  at  which  a  fragment  must  leave  the  shell  in  order  to  hit  the  target  can 
then  be  computed  from  the  following  formula: 

Cos  Ojj  =  sin  u)  cos  t  +  cos  o)  sin  f  cos  $ 

When  the  angle  ©^^  has  been  determined  an  interpolation  in  the  fragmentation 
characteristic  table  is  carried  out  to  obtain  a^,  and  the  count  number 
corresponding  to  the  angle  The  coiont  nmber  is  rounded  off  to  the  nearest 

integer  and  the  associated  mass  distribution  consisting  of  pairs  of  numbers 

+’V» 

m^  and  (i=l,2,w)  (m^  =  representative  mass  of  i  the  weight  group  w  =  nmber 
of  weight  groups,  =  fraction  of  fragments  in  i^^  weight  group)  is  selected. 
Next  the  value  r'  is  computed  from  the  formula: 


6.  Vector  Run 


7.  q  and  m 


8.  Cover 


r' 


/£^JL 

sin  \|f 
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CONCLUSION- 


It  is  hoped  that  the  fairly  detailed  description  of  the  lethal  area 
problem  presented  herein  will  provide  the  reader  with  adequate  background 
not  only  for  the  computation  of  lethal  areas  but  also  to  enable  him  to  consider 
allied  concepts  and  consider  special  cases  of  the  lethal  area  problem.  Since 
virtually  all  of  the  evaluation  models  used  by  the  Weapon  Systems  analyst  are 
based  upon  a  casualty  concept,  i.e.,  expected  fraction  casualties,  number  of 
rounds  or  weight  of  ammunition  required  to  produce  a  given  casualty  level, 
etc.,  an  understanding  of  the  probability  of ■ incapacitation  concept  discussed 
herein  is  mandatory.  Further  many  evaluation  models  utilized  today  are 
approximations  which  utilize  lethal  area  as  such,  giving  some  consideration 
to  weapon  accuracy  and  target  size. 

K.  A.  MYERS 
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